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The Jaynes-Cummings (JC) system, which describes the interaction between a cavity and a two-level atom, is 
one of the most important systems in quantum optics. We obtain analytic solutions for the one- and two-photon 
transport in a waveguide side-coupled to the JC system using input-output formalism in Fock space. With these 
results, we discuss the conditions under which the JC system functions as a photon switch for waveguide photons 
in both the strong and weak coupling regimes 



I. INTRODUCTION 

Controllable interaction between light and matter at the few- 
photon level is one of the main pursuits in quantum information 
processing [1]. Much attention has focused on quantum two- 
level systems in the form of atoms [2] and quantum dots [3- 
5] in the optical frequency range, or Rydberg atoms [6] and 
superconducting Josephson junctions [7-9] in the microwave 
frequency range. All of these different forms of quantum two- 
level systems can be described by the same model; in this paper 
we refer to them simply as "two-level atom". The two-level 
atom is often placed in a cavity, which may enhance or inhibit 
spontaneous emission [10]. In the strong-coupling regime, 
the atom coherently exchanges energy with cavity photons 
[6, 11-15], as observed in vacuum Rabi-splitting [16] of the 
out-coupled photon spectra [4, 17]. 

We consider a system consisting of a waveguide side- 
coupled to a single-mode cavity containing a two-level atom. 
An example of such system is shown in Fig. 1(a) for the case 
where both the waveguide and the cavity are in a photonic 
crystal. This geometry has been demonstrated in several re- 
cent experiments [8, 18-22]. Theoretically, the atom-cavity 
system has been extensively studied with both numerical and 
analytic techniques. Numerically, this system has been simu- 
lated with an inherently stochastic quantum trajectory method 
[23-25], and by directly solving the Master equation with a 
truncated number-state basis describing the photons in the cav- 
ity [18, 26-28]. Analytically, this system can be analyzed 
with a coherent-state input by solving the time-evolution of the 
density matrix [29-31], or by using input-output formalism 
assuming with the weak-excitation approximation [2, 32]. 

We provide a fully quantum mechanical study of one- and 
two-photon waveguide transport through the atom-cavity sys- 
tem. Previously, single-photon transport was studied by solving 
for the real-space representation of the one-excitation eigen- 
state of the system [33, 34], while two-photon transport was 
studied using the field theoretic LSZ formalism [15]; both of 
these solutions are exact. Here, we re-derive the one- and two- 
photon S-matrices using input-output formalism [35] in Fock 
space [36]. our analytic results are in complete agreement with 
Refs. [15, 33]. 
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While the analytic results we present in this paper are known 
in the literature, we believe our re-derivation is of significance. 
First of all, in comparison with the real-space wavefunction, 
and the LSZ field theoretic approach, input-output formalism is 
more widely used in quantum optics. Traditional use of input- 
output formalism, however, has largely focused on solving for 
properties of quantum systems with a coherent or squeezed 
state input. In this regard, it is of value to demonstrate how one 
can use the same formalism to solve for transport properties 
of Fock states in a geometry that is of direct experimental 
importance. Secondly, in the case of two-photon transport, 
input-output formalism in fact results in a much simpler and 
far more transparent derivation of the two-photon scattering 
matrix, compared with the LSZ technique. Finally, we note 
that our Fock-state approach may be generalized for more than 
two photons, and also multiple atoms [37], thus complimenting 
approaches with similar capabilities that are based on coherent- 
state excitations [38-41]. We therefore believe our work here 
can facilitate greater understanding of the transport properties 
of non-classical Fock states in this experimentally important 
geometry. 



In Refs. [28, 42], it was theoretically and experimentally 
demonstrated that the unequal energy level spacing in the JC 
ladder, shown in Fig. 1(b), gives rise to photon-blockade, 
where a single photon blocks the transport of a second, iden- 
tical photon, through a direct-coupled cavity. This effect is 
analogous to the Coulomb blockade effect in condensed matter 
physics, where one electron blocks the transmission of a sec- 
ond electron through sufficiently small semiconductor islands 
due to the single electron charging energy [43]. Using our ana- 
lytic results we show that when the cavity is side-coupled, the 
system can act as a two-photon switch in both the strong and 
weak coupling regimes, and discuss the connection to photon 
blockade in the direct-coupled cavity. 



The structure of the paper is as follows. In Section II the 
Hamiltonian of the system is introduced, and equations of 
motion are derived. In Sec. Ill the single-photon S-matrix is 
calculated, and single-photon transport is discussed. In Sec. 
IV the two-photon S-matrix and the response to a two-photon 
planewave are calculated. In Sec. V we discuss the two-photon 
switch behavior. 
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II. HAMILTONIAN AND EQUATIONS OF MOTION 

A cavity containing a single two-level atom is described by 
the Jaynes-Cummings (JC) hamiltonian [44]: 



Hjc = :^Qcr- + wc"' c + g [: 



c' cr_ + cr.c 



(1) 



where O is the atomic transition frequency, oj is the cavity 
mode frequency, and g is the atom-cavity coupling rate. The 
atom and cavity operators satisfy [cr+, cr_] - cr~ and ^c, c ' j - 1, 
respectively. We side-couple the atom-containing cavity de- 
scribed by Eq. (1) to a waveguide with linear dispersion rela- 
tion and group velocity V;, for right- and left-moving photons 
[Fig. 1(a)], with respective real-space photon creation (anni- 
hilation) operators fl]j(x)[flR(x)] and a|(x)[fl£(x)], satisfying 

[aR(jic), fljj(x')] - [flL(jc), fl|(jc')] - Six - x'). We refer to this 
atom-cavity-waveguide system as the two-mode model. To 
solve the two-mode model, we exploit the spatial inversion 
symmetry of the system by decomposing the Hilbert Space 
into even and odd subspaces with respective photon operators 
adx) = [aR{x) + aii-x)] and Uoix) = [aR{x) - ad-x)], 

satisfying [ae(x), al(jc')j = [ao(x), ao(x')j = 5{x - x'). Both 
even and odd subspaces feature a one-way (chiral) waveguide 
mode, but only the waveguide mode in the even subspace cou- 
ples to the JC system. We refer to the system in the even sub- 
space as the one-mode model. We will solve for the S-matrix 
of the one-mode model, from which the two-mode model S- 
matrix is then straightforwardly constructed (See Appendix A 
in Ref. [36]). The one-mode model hamiltonian is given by: 



i/e — Vg J" dkk aj^a/t + V J" dk |^a|c + c'a^^j -H Hjc, 



(2) 



where V is the waveguide-cavity coupling strength. The 
waveguide photon operator — J dxe''"^ae{x) satisfies 

|^fli,a^j = 5{k - k'), where we omitted the e subscript from 
the A;-space photon operators. Throughout this paper we set 



Heisenberg operator equations for at(f), cr_(f) and c(f) fol- 
low from Eq. (2): 



dt 
dc(0 _ 

"dT ~ 
dcr„(f) 
dt 



-ikakit) - iVc(t), 



I 



iajc(t) - iV I dkokit) - igcr-{t). 



-iQcr_(f) -I- \go-At)c(t). 



(3) 



(4) 



(5) 



In Appendix A we follow the approach of ref. [36], defining 
input and output waveguide photon operators, leading to 



,t(f) = flin(f) - i C(f), 



dc{t) 
dt 



[-10) - ^)c(f) - i V*- flm(0 - cr_(f). 



(6) 

(7) 



where k — 2nV^. Eqs. (5), (6) and (7) are the fundamental 
equations in the input-output formalism for this system. Here 
we will solve these equations for one and two-photon scattering 
matrix (S-matrix). 




FIG. 1. (a) schematic of a two-level atom in a cavity side-coupled to 
a waveguide, (b) the energy spectrum of the IC hamiltonian (fJ = uj) 
in the excitation number manifolds n = 0, 1, and 2. Arrows show 
one and two photon absorption for two-photon switch. The switch is 
observed in reflection in the strong-coupling regime (red arrows), and 
in transmission in the weak-coupling regime (green arrows). 



III. ONE-PHOTON SCATTERING MATRIX 

The single photon S-matrix between two free single-photon 
states with energies k and p may be expressed using input and 
output operators as [36]: 



= {p- I k^) = {0\a,ut{p)ai{km- 



(8) 



Here, \k'^) and \p^) are scattering eigenstates that evolve in the 
interaction picture from free-photon states in the distant past 
and future, respectively [36]. Making use of 



aout(^) = flin(^) - i "V^ c(k). 



(9) 



the Fourier transform of Eq. (6), we rewrite Eq. (8) using 
input-field operators: 

{p-\k^)^5{k-p)-iyrK{mpW)- 

To calculate the matrix element (0 | c{p) \ k^), we derive an 
equation of motion for its time-domain counterpart (0|c(f)l^^) 
using Eq. (7), resulting in 



{Omr) = (-i^^ - ^)<0|c(f)l^^) - i V^<0|ain(f)l^^) 



d 
dt 

- ig{0\cT^{tW) 



(10) 



where <0|ain(f)|/t^) = y^^"'*', and from Eq. (5) (0|cr_(f)|r> 
satisfies 



dt 



{Q\^-{tW) = -iQ(Ok_(0|fc^> - ig{0\c{tW)- (11) 
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Here, we have used the fact that cr,|0) = -|0). Eq. (11) can be 
obtained alternatively using the weak excitation approximation 
by setting cTj — > -1 in Eq. (5). Thus, the weak excitation 
approximation is in fact justified when treating single-photon 
transport. However, in the single photon calculation, the ap- 
proach here can also be used to directly calculate the excitation 
of the atom. 

The solutions to Eqs. (10) and (11) yield: 



(oic(/7)r) - 



= sf6{k-p) 



{0\cT^{p)\e) = 
^sf6{k-p) 



-S(k-p) 



{k-aj + il)(k-n)-g- 



-6{k-p) 



{k-oj-{l){k-Q.)-g^ 

ip-r) = 7 -{ -5{k-p) 

(fe-« + if)(^-n)-^2 

= tk6(k - p) 

where s^'^\k), s^''\k) and fj; are the cavity and atom excita- 
tion amplitudes, and the one-mode transmission coefficient, 
respectively. We note that the poles of all three functions are 



+ g^. In the limit /c — > we have 



Ai -t approaching the one-excitation eigenvalues A^^^ of the JC 
Hamiltonian. 

Having derived the scattering matrix for the one-mode model 
above, we can straightforwardly obtain the transmission and 
reflection amplitudes for the two-mode model: 

R{p\s\kh = mk - p), 
L{p\s\k)R = mk + p). 



where f^. 



(fi + 1), r-k 



(fi - 1) and \k)f.,L 



J iixe'*-*a]j ^(x)|0). Below, we provide a brief discussion of 
the properties of the single-photon transport for the two-mode 
model, while only highlighting those aspects that are relevant 
for the discussions of the two-photon properties. A detailed 
discussion of single-photon transport in this system can be 
found in Ref [33]. 

In Fig. 2, the transmission spectrum for the case of atom 
and cavity on-resonance (a> - Q) is plotted. We plot both the 
weak coupling regime with g < k [Fig. 2(a)], and the strong 
coupling regime with g > k [Fig. 2(b)]. In both regimes, three 
extrema are present in the spectrum. Two transmission min- 
ima where the incident photon is completely reflected occur at 
k - /Ij"^ = ll) + g, the single-excitation energy eigenvalue of 



the JC hamiltonian. At these frequencies, s 



,, the 



cavity and atom excitations have equal amplitudes and either 
equal or opposite phase. . The transmission spectra have a 
maximum at ^ = Q, where the incident photon is fully transmit- 
ted. At this frequency, the cavity excitation amplitude is zero 

(c) 

(s^j^ — 0), while the atomic excitation amplitude is maximal 




FIG. 2. Single-photon transmission IFj^p vs. incoming photon angular 
frequency in the case of tuned atom and cavity (ai = O). (a) strong- 
coupling regime (g > k). (b) weak-coupling regime (g < k). 



s"°l^ — - yfKlg. This is the dipole-induced transparency effect 
pointed out in Ref. [32]. In the strong coupling regime [Fig. 
2(a)], the spectral width of each transmission dip is /c/2, inde- 
pendent of the value of g, while in the weak-coupling regime 
the transmission spectrum displays the electromagnetically- 
induced transparency (EIT) type of feature, with a narrow 
central transmission peak whose width is proportional to g, as 
shown in Fig. 2(b). 



IV. TWO-PHOTON SCATTERING MATRIX 

In the one-mode model, the two-photon S-matrix 

L'e{p\,P2\^ee\k\,k2)ee = (/JiPj l^i^j )' Connecting an incoming 
two-photon state with photon energies k\ and ^2 to an outgo- 
ing two-photon state with photon energies pi and p2, may be 
written as [36]: 

{(i\aout{p2)aout{p\)a]„{k2)a]j,ki)\0) 
= j dp{pi\p^){p^\aUp2)al(k2)al(ki)\0) 

= tp, {pt\aoui{p2)alik2)al{ki)\Q} 

where we have inserted a complete single-photon basis above, 
and made use of the one-photon S-matrix from Sec. III. Sub- 
stituting Eq. (9) we obtain: 

tpi{P'l\am(P2)\klk2} - i ^/K{p^\c{p2)\kik2} 

In order to calculate {p||c(/?2)l^i^2^, we obtain a differential 
equations for its time-domain counterpart {p^\c{t)\kik2) using 
Eq. (4): 



d 
dt 



<K|c(f)l^i^J) = -i(o.-i^)</,||c(f)|fei^J) 
■ i V^<K|fl,„(f)l^ifcJ) - ig<p|k_(f)|fcifc+) (12) 



where {p^\o--(t)\kik2) satisfies 
d 



dt 



<p|k_(0|fcifc2> = -i^{pt\cr-miki) 
■ ig{pt\o-z(t)c(t)\kik^} 



(13) 
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as can be obtained using Eq. (5). We must now solve for the 
matrix element {p^\cr~{t)c{t)\kik2} in Eq. (13), which may be 
rewritten using the identity cr, - 2o-+cr_ - 1 as: 



left-propagating photons following the procedure in Ref. [11]: 

RR{piP2\S\kik2)RR = tkitk2 [(5(^1 - p\)5{ki- p2) + 5(ki - P2)5{k2 - pi)] 



{p\\crAt)cmxkl) = 2{p\\(T4mm(r-{t)c{t)\hkl) 

- {p\m\hkl) (14) 

Examining Eq. (14), we note that the weak-excitation approxi- 
mation (cTj — > -1) certainly does not hold. The first term on 
the right-hand side of Eq. (14), which would be completely 
absent in the weak-excitation limit, is in fact responsible for 
two-photon resonances in the S-matrix. To proceed, we gener- 
ate an equation for the matrix element {Qi\cr_{t)c{t)\k\k^) 



d 

dt 



{0\(T-(t)c(tmk^} = -i + O - i^) <0|£r_(f)c(Ol^i^J) 
- i ^/i{0\cT^{t)aJt)\klk;_} - ig{0\cHt)\kiki} (15) 

where (0|c^(Ol^ifcJ) satisfies 

j^mcHtmk^} = -2i(aj - i'^yO\cHt)\kik^} 

- 2i V^<0|c(f)fl,„(f)|fci^J) - 2ig{0\cTAt)cmik^}. (16) 

In deriving Eq. (16) we have used the identity [c(f), flin(0] - 
(see Appendix A in Ref [37] for a similar proof). We now have 
a closed set of ordinary differential equations for the various 
matrix elements involved. Eqs. (15) and (16) may then be 
solved for {O\o--{t)c{t)\kik^} . Using this solution, we solve 
Eqs. (12) and (13), and obtain the two-photon S-matrix: 



LL{p\P2\S\hk2)RR = fk^rt2 [6{ki + pi)6(k2 + P2) + S(ki + P2)6{k2 + pi)] 
+ \^B5{Ei - Eo) 

RL{p\P2\S\k\k2)RR = tk.fk, [6{ki - pi)6(k2 + P2) + S(ki - p2)6(k2 + pi)] 



+ -B6(E, - Eo) 

We consider an incident two-photon planewave state 
\k\k2)RR, comprised of two right-going photons with individual 
energies k[ and k2, as described by: 



\k\k2)RR = J" dxxdx2Sk^,k2ix\,X2)^i^a\{xi)a\{x2)\0), 

where 5i,,t,(xi,X2) = j^e^i tigazA-z + gai-t2g*2-i;i j a sym- 
metrized two-photon planewave. The resulting outgoing state 
10), calculated in Appendix B, consists of three two-photon 
states 

10) = \4>)rr + \4>}ll + \4>)rl, 

describing two right-moving, two-left moving, and a right and 
left moving photons, respectively, given by: 



1 



\ J, 



{piP2\hK) = hnhi [^(^1 - - Pi) + 5{h - P2)5(k2 - px)] 1^)^^ ^ J dxidx2\TkJk._Su,,k.MuX2) + Hi 



Xl,X2)\^ 



+ B6(Eo - Ei), 



(17) 



V2 



al(xi)a[(x2)\0). 



where 



B 



fkt S k, M ix\ , X2) H- H{xi ,X2)\y. 



(£,-^2,+)(£,-^2,-) 



1 

V2 
\<l>h 



a[{-xi)a[(-X2)\Q}, 



(18) 



(19) 



is the fluorescent term — the source of photon correlation ef- 
fects, and Ei - ki + k2. Eg - p\ + P2 are the total energy of 
the incident and outgoing photons, respectively. The above 
S-matrix is identical to S p^p^k^k^ in Eq. (6) of ref. [15], 
which was obtained using the field-theoretic LSZ method. The 
present derivation is more elementary. Also, in the present 
derivation, the S-matrix' composition in terms of one-photon 
excitation amplitudes and the role of the two-photon poles. 



J dx,dx2^^ [hjk.e^'^-' + hJk^e-^'^-"] + V2i/(xi,X2)J 



a^(ji:i)a[(-X2)|0), 



(20) 



where 



H{xi,X2) = 



F(ki,k2)e' 



lEjXr 



± -I / 1 — I + IS transparent. 



Starting from the S-matrix for the one-mode model in Eq. 
(17), we obtain the scattering amplitudes in the two-mode 
model between planewave states corresponding to right- and 



/(£,-2.i,,-)|Al/2 i(£,.-2ii.+)|x|/2 ] 



{Ei - 2ii,_) {Ei - 2^1,+) 



Above, we have defined the two-photon center-of-mass x^ = 
^^iy^, and spatial separation x = xi - X2 coordinates; we also 
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FIG. 3. Strong coupling regime (g = V5/^). (a) Reflected g<2>(0) 
vs. E/2 in red curve; Transmitted g^^\0) vs. E/2 in green curve, (b) 
Reflected two-photon g'^'Cr) for E/2 = Af ^_^ in red curve; transmitted 
two-photon g*^*(T") for E/2 = Q is green curve. Inset: Single photon 
reflection spectrum Ir^p 



define A,- = 



and 



F(kuk2) = i 



{Ei-A2,^)(Ei-A2,-) 



(21) 



Photon statistics of the transmitted and reflected two-photon 
states is studied through the second-order coherence func- 
tion g^^1(T) = G(2)(t)/|G("(0)|2, where G('>(t) = FF{(f>\al{y + 

T)aF(y)\<p)FF, G<^^\t) ^ff {(l>\al(y)al(y+T)aF(y+T)aF(yM)FF, 
and F - R,L. Here we calculate using the expression 
^'^'(t) = G^'^\t)/G^^\t — > oo), which is a consequence 
of lim^^oo ^'^'(t) = 1. ^^^^(t) is accessible experimentally 
through two-photon coincidence counting. With the help of 
^'^'(t), we now discuss how the JC system can function as a 
two-photon switch. 



V. TWO-PHOTON SWITCH 

Based on the analytic results presented above, we provide 
a discussion of two scenarios where the system can behave 
as a two-photon switch, for energy degenerate photons. Our 
discussion of the strong-coupling regime miiTors that of ref. 
[15]. As a point of departure, here we also discuss a photon- 
switch in the weak-coupling regime, thus providing a unified 
discussion of both regimes. 



A. Strong-coupling regime 

The two-photon switch in the strong-coupling regime is 
a consequence of the anharmonicity of the JC ladder [28]: 
As discussed in Sec. Ill, one incident photon tuned to the 
single-excitation eigenstate of the JC system, k - is fully 
reflected. However, the energy of two such incident photons 
2k - Ej - 2/lj°^, does not coincide with a two-excitation eigen- 
state of the JC system, as illustrated by red arrows in Fig. 1(b). 
Consequently, the two photons are prevented from simultane- 
ously entering the cavity. Since reflection in the side-coupled 
cavity considered here arises entirely from cavity field decay, 
sub-poissonian statistics g'^^(O) << 1 is observed for the two 
reflected photons, as shown in Fig. 3(a). Additionally, g^^Xr) 
of the two reflected photons exhibits strong anti-bunching and 
sub-poissonian statistics for all time intervals r, as shown in 
Fig. 3(b). Here, a sufficiently large coupling rate g is needed 
to ensure that the two-photon energy does not overlap with the 
two-excitation energy eigenstates, which are broadened due to 
waveguide coupling. We note that the photon switch behavior 
in reflection described above in the side-coupled cavity is iden- 
tical to the photon-blockade reported in transmission through 
a direct-coupled cavity [28, 42]. 

In addition to the photon switching behavior in reflection 
for Ei/2 - /Ij"^, the g^^'(O) spectrum also exhibits other in- 
teresting features. In reflection, g*^*(0) peaks at Ej/2 - Q, 
exhibiting super-poissonian statistics [Fig. 1(a)]. This be- 
havior is attributed to the vanishing single-photon reflection at 
^ = Q, since §'^*(t) is normalized with respect to G*''(0), which 
vanishes when single-photon reflection vanishes. We note 
that in the direct-coupled cavity system, the super-poissonian 
statistics shown here in two-photon reflection at Ei/2 — Q 
occurs in transmission instead [42, 45], and is named 'photon- 
induced tunneling' [42]. In transmission, g^^\Q) peaks at 
Ei/2 - /Ij"^, displaying super-poissonian statistics [Fig. 1(a)], 
which is attributed to the vanishing single-photon transmission 
at k\ - Finally, in transmission, g^^\Qi) ~ 1 for all time 
intervals t at Ei/2 - Q.. At this frequency, the atom excitation 
scales as 1 /g, and is therefore weakly excited under strong- 
coupling, resulting in poissonian statistics for the transmitted 
light [Fig. 1(a)]. 



B. Weak-coupling regime 

A two-photon switch in the weak-coupling regime is also 
possible in this system, without requiring a nonlinear cavity 
[46]. Here, the photon switch is created due to the intrinsic 
non-linearity of the atom [47, 48]. As discussed in Sec. Ill, 
one incident photon tuned io k — Q., is completely transmit- 
ted. However, when two such photons, each with an energy 
Ei - 2n, are incident, sub-poissonian statistics g^'^\0) « 1 
is observed for the two transmitted photons, as shown in Fig. 
4(a). g'^^\T) of the two transmitted photons is anti-bunched, 
and displays sub-poissonian statistics for all time intervals t, 
shown in Fig. 4(b). Thus, the presence of the atom in the 
cavity results in complete single-photon transmission (dipole- 
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induced transparency), but the inability of the atom to absorb 
more than one photon at a time leads to anti-bunching and sub- 
poissonian statistics in two-photon transmission. We note that 
a similar atom-induced photon switch occurs when the atom is 
directly coupled to a waveguide [11]. There, a single-photon 
is completely reflected, and anti-bunching is observed in two- 
photon reflection g^^\0) - 0, while bunching is observed in 
two-photon transmission. 

Moreover, the two-photon reflection at E,/! = Q displays 
super-poissonian statistics as a result of the vanishing single- 
photon reflection at A; = Q. Finally, as is the case of strong- 
coupling, super-poissonian statistics is observed in two-photon 
transmission at Ej/l - Af]. 



VI. FINAL REMARK AND CONCLUSION 

As a final remark, we note that in practice, in addition to cou- 
pling to the cavity which is in turn coupled to the waveguide, 
the atom also couples to non guided modes, leading to loss. 
This loss may be accounted for by making the replacement 
Q — > Q - iy/2 where y is the coupling rate into non-guided 
modes [49]. For the case of the D2 transition of a Caesium 
atom placed in a Fabry-Perot cavity, this loss can be on the 
order of the cavity linewidth [28]. In contrast, for the case of a 
self-assembled InAs quantum dot in a photonic crystal cavity, 
the atom loss is roughly two orders of magnitude smaller than 
the cavity linewidth [42], and therefore does not lead to a qual- 
itative difference in the transport properties with respect to the 
non-lossy case. Here, we therefore ignore the effects of atomic 
loss. 

In conclusion, we have solved for the one- and two-photon 
S-matrices in a waveguide side-coupled to a cavity containing 
a two-level atom. Our solution is based on input-output formal- 
ism [35], in Fock-space [36], and is fully quantum-mechanical 
and deterministic. We have discussed the features of one- 
photon transport, and the photon-switch effect in both the 
weak- and strong-coupling regime, as seen in the transport of a 
two-photon planewave state. Using the S-matrices presented 
in this paper, the response to one and two incident waveg- 
uide photons with arbitrary spectra and pulse shapes may be 
straightforwardly obtained. 



Appendix A: Derivation of input-output relations for waveguide 
photon operators 

Following the approach of ref. [36], we integrate Eq. (3) 
from time fo to obtain 



(a) 



iV f dt'c(t')e-''''^'-''> (Al) 

J to 



We define k = InV^, the waveguide-cavity coupling rate. Ad- 
ditionally, we define the operator 0(f) = J dkaic(t), the 

input-field operator a,„(f) = ^ J t/^e"'**'"'°'a<:(fo), and sum 




FIG. 4. Weak coupling regime (g = -^). (a) Reflected g*^'(0) vs. 
E/2 in red curve; Transmitted g*^'(0) vs. E/2 in green curve, (b) 
Reflected two-photon g'''(T) for E/2 = Af ], in red curve; transmitted 
two-photon g'-^\T) for E/2 = Clis green curve. Inset: Single photon 
transmission. 



Eq. (Al) over k, leading to: 



<D(f) = fli„(0-iA/^c(f) 



(A2) 



We also define an output field operator aout(t) = 



J dke '''a^(fi), integrate Eq. (3) from time fj oo. 



and sum over k, leading to: 



(A3) 



^(t) = Oontit) + i J -C(t) 

Equating Eqs. (A2) and (A3) we obtain: 

flout(f) = flin(f) - i yfKc(t). 

Plugging Eq. (A2) into into Eq. (4), we obtain: 

dc / r 

— = - I2 jc(0 -i^KOmit) - igo--(t). 

We have thus obtained the input-output formalism [Eqs. (6), 
(7)] for this system. 



Appendix B: Two-photon planewave transport 

Consider an incident two-photon planewave state with right- 
moving photons: 

\k\,k2)RR = J" dxidx2Sk^,k2i^\,X2)-^a\(xx)a\{x2)\Q\^\) 
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where St^M^ (xi , X2) = ^ [e*' " e'*^"^^ + e""'''-e'''- '" ] . We de- noting that eiplS^ \k)e = h6(k - p) we have 
compose the state in Eq. (Bl) into even and odd subspaces by 



2 V2 



\xi,X2}rr = ^fljj(xi)a^(x2)|0) 



^ ^ fal(xi) + fl^(xi)] [cl(x2) + 4(X2)1 10) 



1 1 

1 



al(x,)al,{x2m. 



-al{xi)al{x2)\Q) + ^fll(xi)fl^(x2)|0) + -^al{x^)al{x2W'^''^^^' ^" *^ subspace we have: 
2V2 2V2 2V2 

1 



+ -^al{x{}al{x2)\0) 
2 V2 

= \\Xl,X2)ee + \\x\,X2)oo + -^\x\,X2)eo + -^\xi,X2)o 

2 2 2V2 2V2 



2 V2 



SJ|^l,^2>o 



— dxidx2— — 



2V2 

al(xi)al(x2m 



It follows that 

2 2 2 V2 



1 



2 V2 



\kl,k2}oe. 



where \ki,k2}ee — J dxidx2S k,,k2(xuX2)-^al(xi)al(x2)\0) and 

(e.g) \ki,k2}eo = / dxidx2S k,,k2(xu X2)al,(xi)a],(x2)\0}. We can 
now apply the scattering operator: 

S\kuk2)RR = l;See\kl,k2)ee + ^Sool^i , ^2)00 + ^-pS^l^i , ^2>co 

2 2 2 V2 



+ ]=Soe\ki,k2)o 

2 V2 



(B2) 



where S|^i, ^2)«r is the out-state (scattered state). Below, we 
work through each one of the terms in Eq. (B2). 



1. eo and oe subspaces 

In the eo subspace the scattering operator is S^o - S^So = 
Selo, where !„ is the identity operator in the odd subspace. We 
may then rewrite the third term in Eq. (B2) 



1 



Selo\kuk2)eo = 



2 V2 

= — f dxidx2—^\e"""e'''''' +e""''e"''"]x 
2 V2 J V227r ^ ^ 

S,al(xi)al(x2m. 

Inserting a one-photon resolution of the identity twice, we 
rewrite as: 

= — r dxidx2—^\e""''e'''^''- +e"""^e"'^"]x 
2 V2 J V22;r ^ ^ 

J dk J dp\p)e e{p\'&e\k)e e{k\x\)e\X2)o 



2. 00 subspace 

In the 00 subspace the scattering operator is Son = SqSo = 
IqIo where is the identity operator in the odd subspace. 
Applying the scattering operator we have: 



^SoSo\ki,k2)oo = ^1^1,^2)0 



3. ee subspace 

In the ee subspace we have 

1 1 r°° CP' 

-See\kl,k2}ee ^ dp\ dp2\p\, P2)ee ee{P\, P2\^ee\k\,k2)ee, 



where we have inserted a two-photon resolution of the identity. 
We use the result in Eq. (17) 



-See\kl,k2)ee ' 2 j dp\dp2\tkitk2 



6(ki - p\)6{k2 - P2) 



+ 6(ki -P2)6(k2-pi) 



\PuP2)e 



+ sf>sf^F(ki , ^2)5(^1 +k2-pi-p2)\x 



(B3) 



where \pi,p2}ee - ^ dxidx2e'^"''' cos ls.oX\xi,X2)ee, and 
F(k\,k2) was defined in Eq. (21). The first term in Eq. (B3) — 
the uncorrelated transport term — is equal to \tkjk2\ki,k2)ee- 
The second term is equal to 

^ r 

— F{kuk2) I c/xiliX2|xi,X2)ce 

dEj^-''6{Ei - £„) d^„ cos \^oX\sf^sf^. 

DO xj — OO 
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Integrating over £„ and exploiting the invariance of the inte- Integrating over Ao using the residue theorem: 
grand under the inversion A,, — » -Ao: 



F(ki,k2) J dxidx2e'^'''''X 



2V2 

j(E,/2-At,_)\x\_ 1 1 



_^ gi(£,/2-i,.+)|.t| 



1 1 



Finally, we combine the uncorrelated and correlated terms 
and get 

-See\kl,k2)ee - ^ '/ti 'fc 1^1 ' ^2)cc 



Or: ig2 



+ 



• 2 n 

2 V2 J 



2 V2 

g;(£,72-^,,_)|x| 1 1 



«^ — cx 



^ gi(£,72-^i,+ )|x| 

-F(kuk2) ^ dxidX2e''^''''\xi,X2)ee \xi,X2)ee- 



- ^1,-) {Ei - 2ii,_) (£,■ - - 

1 1 



^2^2^ ^ ('ti,- - ^^1,+) - 2/li_+) (£■/ - - /li _) 
Stt 

i^Ape'^"-^ ^ We then add together the results of sections B 1 , B 2 and B 3, 

(A„ + EJ2 - ^,,_) (A„ + £,/2 - A,,,) ^^^^^^-^^ t ^ 4 = to 

V2 " V2 

express the out-state in terms of left and right moving photons, 



1 



(Ao - £,/2 + (Ao - £,/2 + ^1,+) as presented in Eqs. (18), (19) and Eq. (20) 
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